clc,clear
close all
F =[];
IndexF = 1;
parfor RotatetTheta = 0:1:30
    F = [F FCaseFan(RotatetTheta)];
end
% load("数值积分")
% plot([0:1:30],F(1);
% xlabel("\theta/deg");
% ylabel("Tz/(N*m)");